• Journal of Mathematical Physics 51, 113504-10 (2010) • 



A CLOSED FORMULA FOR THE BARRIER TRANSMISSION COEFFICIENT 
IN QUATERNIONIC QUANTUM MECHANICS 



Stcfano De Leo 1 , Gisele Ducati 2 , Vinicius Lconardi 1 , and Kenia Pereira 1 . 

1 Department of Applied Mathematics, State University of Campinas, 
SP 13083-970, Campinas, Brazil, 

deleo Qime.unicamp.br. 

2 CMCC, Federal University of ABC, 
SP 09210-170, Santo Andre, Brazil, 
ducati&ufabc. edu. br. 

Submitted: May, 2010. Revised: Sepetember, 2010. 

Abstract. In this paper, we analyze, by using a matrix approach, the dynamics of a non- 
relativistic particle in presence of a quaternionic potential barrier. The matrix method 
used to solve the quaternionic Schrdinger equation allows to obtain a closed formula 
for the transmission coefficient. Up to now, in quaternionic quantum mechanics, almost 
every discussion on the dynamics of non-relativistic particle was motived by or evolved 
from numerical studies. A closed formula for the transmission coefficient stimulates an 
analysis of qualitative differences between complex and quaternionic quantum mechanics, 
and, by using the stationary phase method, gives the possibility to discuss transmission 
times. 
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I. INTRODUCTION 



In the quaternionic formulation of non-relativistic quantum mechanics, the dynamics of a particle 
without spin subject to the influence of the anti-hermitian scalar potential, 



iV 1 (r,t)+jV 2 (r,t) + kV 3 (r,t) , 



is described by 



hd t ^(r,t) = 



h 2 

2m 



iV 1 (r,t)-jV 2 (r,t)-kV 3 (r,t) 



*(»■,*) 



(1) 



Eq.(l) is known as the quaternionic Schrodingcr equation. For a detailed discussion on foundations of 
quaternionic quantum mechanics, we refer the reader to Adler's book [1]. Numerical investigations on 
the observability of quaternionic deviations from the standard (complex) quantum mechanics [2-10] 
have recently stimulated the study of new mathematical tools in solving quaternionic differential 
equations [11-16]. 

In complex quantum mechanics [17,18], the rapid spatial variations of a square potential intro- 
duce purely quantum effects in the motion of the particle. Before beginning our analytic study of 
the quaternionic potential barrier and analyze analogies and differences between the complex and 
quaternionic dynamics, we introduce some important mathematical properties of the quaternionic 
Schrodingcr equation in the presence of time independent potentials, 



hd t V(r,t) 



i V -iV^-jV^-kVsir) 
Am 



(2) 
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Due to the fact that our analysis is done for a time independent potential, we can apply the method 
of separation of variables. So, we introduce 



Hf(r,t) = $(r) exp[-iEt] 



(3) 



The position on the right of the complex exponential is important to factorize the time dependent 
function. Then, Eq.(2) reduces to the following quatcrnionic (right complex linear) ordinary differential 
equation [12] 



E${r)i 



i — V-iVAr)- j Voir) - k V 3 (r) 
2m w 



$(r) . 



(4) 



In the case of a one-dimensional potential, V(r) = V(x), and for a particle moving in the x-direction, 
$(r) = $(x), the previous equation becomes 



E$(x)i 



i — d xx -i V 1 {x) - j V 2 (x) - k V 3 (x) 



<f>(x) 



(5) 



In the case of a quaternionic barrier, 

V(x) = { zone I : (x < 0) , zone II : V (0 < x < L) , zone III : (x > L) } , 

the plane wave solution is obtained by solving a second order differential equations with (left) constant 
quaternionic coefficients [12,16]. To shorten our notation, it is convenient to rewrite the quaternionic 
differential equation (5) in terms of adimcnsional quantities. In the potential region, < £ < A, the 
differential equation to be solved is 



* <(0 = * V c ® u (0 + j V q e~ i6 $ n (C) - e 2 i 



(6) 



where 



= ^EjV, , V = yjv 2 +V 2 2 + V 2 , V c = V/V , V q = yJv 2 +V 3 2 /V , arctanfl = V 3 /V_ 



and 



£ = ^2mV /h 2 x , A = ^2mV /h 2 L . 



The general solution contains four complex coefficients (A, A, B, and B) to be determined by the 
boundary conditions, 

*n(0 = (1+J7) {cxp[a_£] A + cxp[-a_£] B} + ([3 + j) {exp[a + £] I + cxp[-a + £] s} , (7) 

for a detailed derivation of this solution see ref. [16]. The potential and energy inputs appear in the 
following complex quantities 



V c ± 



V 2 







V q e 



V q e 



e J - V„ 



e 4 - V„ 



(8) 



Note that for E > ^JV 2 + V 2 , we have 7 = 0*. 

In the free potential regions, zone I (£ < 0) and zone III (£ > A), the solutions are given by 

$:(£) = cxp [ i e£] + exp [ — ie£] R+J cxp [e£] R , 
$ m (0 = exp[ie£] T + j exp[-e£] T . 



(9) 



It is important to observe that starting from theses solutions we can find the relation between reflec- 
tion and transmission coefficients without the necessity to know their explicit formulas. In fact, for 
stationary states the density probability is independent from the time and, consequently, the current 
density <&t(x) i <&'(x) + h.c. is a constant in x. This implies, for example, 



$t(0) z $1(0) + h.c. = $t n (A) 1 $; n (A) + h.c. 
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From this relation, we get the well known norm conservation formula, i.e. 

1-\R\ 2 = \T\\ (10) 

Having introduced this preliminary discussion on the quaternionic Schrodinger equation, we can now 
return to our task of finding a closed formula for the transmission coefficient. Once obtained the 
transmission coefficient it suffices to use Eq.(10) to find the norm of the reflected wave. 



II. CONTINUITY CONSTRAINTS AND MATRIX EQUATION 

To find a closed formula for the transmission coefficient, we impose the continuity of <!?(£) and its 
derivative, at the points where the potential is discontinuous, i.e. 

* x (0) = MO) , 3(0) = Mo) , *m(A) = MA), M(A) = MA). 

Using these constraints and separating the complex from the pure quaternionic part, we find 

1 + R = A + B + 0(A + B) 

A-B -/3(A-B) 

B , 



(l-R) 
R 



P(A 
o+ 

a_ 

1 {A + B) + A 



^R 

o_ 

ieX 



j(A — B) + — (A — B) 



Te 
i— Te 1 



Ae 



Te 

e ~ 
— Te 



- c A 



+ Be 

e 

a_A . 



a_ A 



j3(Ae 



CK_|_ A 



Ae a ~ x -Bc- a - x 



Be; 

-A 



■ a+ An 



— /3(ie a+A -Be- a+A ) 



7(ie 
1 (Ae a ~ x 



Be 
Be 



-a_A 



— a_ A 



) + Ae 



a+ A 



Be 



- a_L A 



) + — (AS 

a_ 



Be 



— CK-L A 



(11) 



The procedure for calculating reflection and transmission coefficients by using the transfer matrix 
method is a well-known technique both in quantum mechanics [17, 18] and optics [19,20]. It is based 
on the fact that, according to Schrodinger or Maxwell equations, there are simple continuity conditions 
for the field across boundaries from one medium to the next. If the field is known at the beginning of a 
layer, the field at the end of the layer can be derived from a simple matrix operation. The final step of 
the method involves converting the system matrix back into reflection and transmission coefficients. 
From Eq.(ll), after straightforward algebraic manipulations, we find 



where 





( 


1 


F = 


i 


a - 






{ 






R 1 

R 
\RJ 



F 




and 



diag{ e 



-a_A 



G 



_a_ A 



i 

■ a+A 



e 4eA ^ 

gieA 










e -eX 






e~ eX J 






1 




P\ 


-1 2± 




7 1 

7 it 


1 

a- / 



(12) 



- e 



a i A 



} 



From this matrix system, we can extract the transmission coefficient T. The closed formula is not 
simple and we prefer to give it in terms of the matrix elements of M, 



T = 2 exp[- ie\]/V , 



(13) 
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where T> is given by 



. e 2 M 12 -a 2 _M 21 



e o_ 



M 13 + iM 24 - 



e 2 M li + ia 2 _M 23 



e a_ 



M 31 - iM 42 + 



i e 2 M 3 



A'L 



M 44 + M 3[ 



e 2 M 34 + a 2 _M 4 , 



The explicit formulas of the elements of the matrix M are 



(1 


— zw) 


Af u = 


+ cosh(a_ A) — 


,87 cosh(a + A) , 


(1 


— zw) 


M 12 = 


— sinh(a_A) + 


a_ + /?7 sinh(a + A) 


(1 


— zw) 


M 13 = 


— /3 [cosh(a_ A) 


— cosh(a + A) ] , 


(1 


— zw) 


M 14 = 


+/3 [sinh(a_A) 


— a_ + sinh(a + A) 


(1 


— zw) 


M 21 = 


— sinh(a_ A) + 


a + _ /?7 sinh(a + A) 


(1 


— zw) 


M 22 = 


+ cosh(a_ A) — 


/37cosh(a + A) , 


(1 


— zw) 


M 23 = 


+f3 [sinh(a_A) 


— a + _ sinh(a + A) 


(1 


— zw) 


M 24 = 


— [cosh(a_ A) 


— cosh(a + A) ] , 


(1 


— zw) 


M 31 = 


+7 [cosh(a_ A) 


— cosh(a + A) ] , 


(1 


— zw) 


M 32 = 


—7 [sinh(a_A) 


— a_ + sinh(a + A) 


(1 


— zw) 


M 33 = 


— f3-f cosh(a_A) 


+ cosh(a + A) , 


(1 


— zw) 


M 34 = 


+(3j sinh(a_A) 


— a_ + sinh(a + A) 


(1 


— zw) 


M 41 = 


—7 [sinh(a_A) 


— a + _ sinh(a + A) 


(1 


— zw) 


M 42 = 


+7 [cosh(a_ A) 


— cosh(a + A) ] , 


(1 


— zw) 


Af 43 = 


+(3j sinh(a_A) 


— a + _ sinh(a + A) 


(1 


— zw) 


M 44 = 


— /?7 cosh(a_ A) 


+ cosh(a + A) , 


Q 


-+ = a 


-/<*+• 







(14) 



where a + _ = a + /a_ 

A first important result is the evidence that the transmission coefficient does not depend on 6. 
This not trivial result is due to the fact that the closed formula obtained for T contains the terms 



M„ 12 



and AL 



33,34,43,44 



which arc independent of 9, their formulas contain the quantity ,87 which 



is independent on 6. and the terms proportional to /3, i.e. Af 13 . 14 . 23 . 24 , which are multiplied by terms 
proportional to 7, i.e. Af 31i32i41j42 . The invariance of the transmission coefficient in the plane (V 2 , V 3 ), 
seen in the numerical studies appeared in littcraturc [7], is now proved to be a property of quaternionic 
quantum mechanics. 

The results of complex quantum mechanics can be obtained by taking the limits (V c , V q ) —> (1,0). 
Observing that a ± — > y/1 ± e 2 and (z, w) — > (0, 0), we find 



M 12 
M 13 
M 14 
M 21 = 
Af 22 = 
M 23 = 
M 24 = 
M 31 = 
M M = 



+ cosh( yl 
— sinh(v / l _ 
, 
, 



j 2 A) 
~?A) 



sinh(V± - e 2 A) 
cosh(Vl - e 2 A) 



, 
, 
. 

, 

M 33 = +cosh(vTTe 2 A) , 

Af 34 = -y/(l-e')/(l + e') sinh(vTTFA) , 

M 41 = , 

M 42 = , 

Af 43 = -y/(l + e')/(l =F) sinh(vTTFA) , 
Af 44 = +cosh(Vl + e 2 A) . 



In this limit, the expression for T> reduces to 



V c = 2 cosh(Vl - e 2 A) + i 



1 - 2e 2 



(15) 



sin 



h(v / T _ e 2 A) , 
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and the transmission coefficient becomes [17, 18] 

7 



T 



— ie\ 



cosh(Vl - e 2 A) +i 



1 - 2e 2 



2ey/T^e' 



sinh(Vl - e 2 A) 



In the case of diffusion, e > 1, the transmission coefficient is often given in terms of cosine and sines 
functions, i.e. 

1 - 2e 2 



Tr = e 



— ie\ 



cos(Ve 2 - 1 A) +i 



2ey/e 2 ^T 



sin(Ve 2 - 1 A) 



(16) 



III. RESONANCE'S PHENOMENA 

In this section, we analyze when possible analytically and when not numerically the phenomenon of 
resonances. For the complex case, the transmission probability 



1 



sin 2 (V^TA) 



(17) 



4e 2 (e 2 - 1) 

shows that the phenomenon of resonances, \T C \ = 1, happens when the energy /barrier width condition 

x/e 2 n -l \ n =n% (18) 
is satisfied [17,18]. The transmission probability oscillates between one and 

1 



1 + 



4?„ 2 (?„ 2 -l) 



(19) 



obtained for \J e 2 — 1 X n = (2 n+1) it / 2. This minimum tends to one for increasing incoming energies. 

• Fixing the width of the potential barrier, A„ = A , and varying the energy of the incoming particle, 
we can find the energy values, (e„,e„ + Ae„), for which the transmission probability reaches two 
consecutive maxima, and the energy value, e„ + A<f„ of the minimum between such maxima, 



e 2 n -l = n 2 ir 2 1 A" , 
(e„ + Ae n f - 1 = (n + 1/2) V /A* 
(e„ + Ae„) 2 - 1 = (n + l)\ 2 / \ 2 



Ae r 



y /l + (n + 1/2) 2 ^ 2 / Ag - y/ 1 + n tt 2 / A 2 
Ae„ = y/l + {n + I) 2 tt 2 I X 2 - y/l + n 2 n 2 / X 2 . 



For complex potential, (V c , V q ) = (1, 0), and for a barrier width of 37r, we find 



e 1 = VlO/3 - 1.054 , Ae x = (V13- \/I6)/3 - 0.148 , 
e 2 = \/l3/ 3 - 1.202 , Ae 2 = y/2 - V13/ 3 - 0.212 , 
e 3 = \/2~ 1-414 . 

Fixing the energy of the incoming particle, e„ = e , and varying the barrier width, we find 



A„ = n 7r / yje\ - 1 



A„ + AA = (n + 1/2) Tr/ Ve 2 - 1 
A„ + AA= (n + 1) tt / V e 2 - 1 



AA = 7r/2 v^= T, 
AA = tt/ . 



For complex potential, (V c , V q ) = (1, 0), and for an incoming energy e = y/2, we find 

A 1= 2 7T , AA!=7T, 
A 2 = 3 7T , AA 2 = 7T , 
A 3 = 4 7T . 
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In Fig. 1, where the barrier width is set to 3tt, we see the phenomenon of resonance for complex, 
mixed and pure quaternionic potentials varying the energy of the incoming particle. It is interesting 
to observe that e„ and Ae n decrease when the quaternionic part of the potential tends to one, 



(Vc,V q ) 




e 2 


A 6l 


e 3 


Ae 2 


(1,0) 1.054 


1.202 


0.148 


1.414 


0.212 


(V3/2,l/2) 1.049 


1.188 


0.139 


1.394 


0.206 


(1/V2,1/V2) 1.043 


1.170 


0.127 


1.369 


0.199 


(l/2,V3/2) 1.034 


1.145 


0.111 


1.334 


0.189 


(0,1) 


1.011 1.077 


0.066 


1.246 


0.169 



In Fig. 2, where the incoming energy is set to v2, we see the phenomenon of resonance for complex, 
mixed and pure quaternionic potentials varying the barrier width. It is interesting to observe that AA 
is constant for a fixed potential and decreases when the quaternionic part of the potential tends to 
one, 



(V c ,V q ) 


\M 


A 2 H 






AA 2 [ttJ 


(1,0) 2 


3 


1 


4 


1 


(V3/2, 1/2) 


1.949 


2.915 


0.966 


3.881 


0.966 


(1/V2, 1/V2) 


1.890 


2.817 


0.927 


3.744 


0.927 


(1/2, V3/2) 


1.819 


2.695 


0.876 


3.571 


0.876 


(0,1) 


1.718 


2.478 


0.760 


3.238 


0.760 



IV. THE LIMIT CASE BETWEEN DIFFUSION AND TUNNELING 

In this section, we present the interesting case of e = 1, i.e. the limit situation between diffusion and 
tunneling. We start by analyzing this limit for a complex potential (V c , V q ) = (1, 0). In the free region, 
the solution of the Schrodingcr equation, for e = 1, becomes 

*' I (0=eacp[ie] + exp[-ie] K+3 cx P [£] R' c , 
<&• m (0 = exp [ it } T' + j exp [ - £ ] f' c . 

Eq.(6) reads 

* % *• n (0 = * Ku (£) - n (0 i, < e < A . (20) 
Separating the complex from the pure quaternionic part in the solution, i.e. 

we find two complex differential equations which lead to 

€,n(£) = C' c exp[V2t] + D' c cxp[- v^C] , (21) 

observe that only two coefficients appear in the complex solution ip* n (£). Imposing the continuity 
conditions, we find 

R' c = - i A / (2 - i A) , 

T' = 2 exp[ — i A ] / (2 — i A) , (22) 

and as expected R' = T* = reproducing the standard complex dynamics. In the case of pure 
quaternionic potentials, the free solutions are 

=exp[i£]+exp[-i£] R'+j exp[t] R' q , 
^,m(£)=exp[i£] r;+jexp[-e] f' q , 
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and the equation to be solved is 

* % s; n (0 = j e- ie *; n (o - *; n (o %, o < e < a . 

Separating the complex from the pure quaternionic part, 

we find, as solutions of the two complex differential equations coming from Eq.(23), 



(23) 
(24) 



(25) 



By imposing the continuity conditions, we obtain the reflection and transmission amplitudes for pure 
quaternionic potentials in the limit case between diffusion and tunneling, 



R' = - iX (6 + 4A + A")/[24 + 24(l-i)A- 18 z A — 4 (1 + i) A — A ] , 

T' = 2 expHA] (12 + 12A + 6A 2 + A'') / [24 + 24 (1 - i) A - 18 i A* - 4 (1 + i) A* - A 4 



(26) 



The expression obtained for T* is in perfect agreement with the limit for e — > 1 of the coefficient T 
given in Eq.(13). 

The qualitative differences between complex and pure quaternionic potentials can be, for example, 
seen comparing the reflection and transmission coefficients in the case of thin (A <C 1) and large 
(A 1) barriers. The computation taking into account fourth order corrections gives 



A « 1 



A> 1 



R'\ ~ A/2- A 3 /16 + 0[A 5 



R' q \ 



R'd 



R' q \ 

T' q \ 



1 - A"/8 + 3A /128 + 0[A' ] , 

A 2 /4- A 3 /12 + 0[A 5 ] , 
1- A 4 /32 + 0[A 5 ] , 

1-2/A 2 +6/A 4 +0[1/A 5 ] , 
2/A-4/A 3 +0[1/A 5 ] , 

1 - 2/A 2 - 8/A 3 + 6/A 4 + 0[ 1/A 5 ] , 
2/A + 4/A 2 - 8/A 3 - 8/A 4 + 0[ l/A* ] 



V. CONCLUSIONS 



With the development of a consistent theory of quaternionic differential equations [11-16], it is now 
possible to be more profound in discussing quaternionic quantum mechanics. The problem of diffusion 
by a quaternionic potential, limited up to now to numerical analysis [2-10], was solved by a matrix 
approach and leaded to a closed formula for the transmission amplitude. 

As in the case of standard quantum mechanics, for incoming particle of energy e > 1, quaternionic 
barriers present the phenomenon of resonances. A comparison between the complex, V c = 1, and the 
pure quaternionic, V q = 1, cases show that, for pure quaternionic potentials, the minimum value of 
transmission increases whereas the oscillation interval decreases. 

This paper could motivate the study of quaternionic quantum mechanics by using the wave packet 
formalism. An interesting application is the study of tunneling times [21,22]. The Hartman effect [23] 
surely represents one of the most intriguing challenges recently appeared in literature [24]. Quater- 
nionic deviations from the standard phase, which is fundamental in the calculation of tunneling times, 
can be now investigated by using the phase of the transmission coefficient given in this paper. 

The analogy between quantum mechanics and the propagation of light through stratified me- 
dia [25] suggests to analyze in details the limit case e = 1. In this case the momentum distributions 
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should be centered in E a = V and phenomena of diffusion and tunneling must be treated together. 
The outgoing wave packets should move with different velocity and should give life to new interesting 
phenomena of interference. 

Acknowledgements. One of the authors (SdL) wish to thank the Department of Physics, University 
of Salcnto (Lecce, Italy), where the revised version of the paper was prepared, for the invitation and 
the hospitality. SdL also thanks the FAPESP (Brazil) for financial support by the grant n. 10/02216-2. 



References 

1. S. L. Adler, Quaternionic quantum mechanics and quantum fields (New York, Oxford University Press, 
1995). 

2. A. Peres, "Proposed test for complex versus quaternion quantum theory", Phys. Rev. Lett. 42, 683-686 
(1979). 

3. H. Kaiser, E. A. George and S. A. Werner, "Neutron interferometric search for quaternions in quantum 
mechanics", Phys. Rev. A 29, 2276-2279 (1984). 

4. A. G. Klein, "Schrodinger inviolate: neutron optical searches for violations of quantum mechanics", Phys- 
ica B 151, 44-49 (1988). 

5. A. J. Davies and B. H. McKellar, "Non-relativistic quaternionic quantum mechanics", Phys. Rev. A 40, 
4209-4214 (1989). 

6. A. J. Davies and B. H. McKellar, "Observability of quaternionic quantum mechanics", Phys. Rev. A 46, 
3671-3675 (1992). 

7. S. De Leo, G. Ducati and C. Nishi, "Quaternionic potential in non-relativistic quantum mechanics, J. 
Phys. A 35, 5411-5426 (2002). 

8. S. De Leo and G. C. Ducati, "Quaternionic bound states, J. Phys. A 38, 3443-3454 (2005). 

9. S. De Leo and G. C. Ducati, "Quaternionic diffusion by a potential step, J. Math. Phys. 47, 082106-15 
(2006). 

10. S. De Leo and G. C. Ducati, "Quaternionic wave packets, J. Math. Phys. 48, 052111-10 (2007). 

11. S. De Leo and G. Scolarici, "Right eigenvalue equation in quaternionic quantum mechanics", J. Phys. A 
33, 2971-2995 (2000). 

12. S, De Leo and G. C. Ducati , "Quaternionic differential operators, J. Math. Phys. 42, 2236-2265 (2001). 

13. S. De Leo, G. Scolarici and L. Solombrino, "Quaternionic eigenvalue problem", J. Math. Phys. 43, 5812- 
2995 (2002). 

14. S, De Leo and G. C. Ducati , "Solving simple quaternionic differential equations, J. Math. Phys. 44, 
2224-2233 (2003). 

15. S, De Leo and G. C. Ducati , "Real linear quaternionic differential operators, Comp. Math. Appl. 48, 
1893-1903 (2004). 

16. S, De Leo and G. C. Ducati , "Analytic plane wave solution for the quaternionic potential step", J. Math. 
Phys. 47, 082106-15 (2006). 

17. C. Cohen- Tannoudji, B. Diu and F. Laloe, Quantum mechanics (New York, John Wiley & Sons, 1977). 

18. E. Merzbacher, Quantum mechanics (New York, John Wiley & Sons, 1970). 

19. O. S. Heavens, Optical Properties of Thin Films (Butterworth, London, 1955). 

20. M. Born and E. Wolf, Principles of optics: electromagnetic theory of propagation, interference and diffrac- 
tion of light (Oxford, Pergamon Press, 1964). 

21. E. H. Hauge and J. A. Stovneng, "Tunneling times: a critical review", Rev. Mod. Phys. 61, 917-936 
(1989). 

22. V. S. Olkhovsky, E. Recami, and J. Jakiel, "Unified time analysis of photon and particle tunneling" , 
Phys. Rep. 398, 133-178 (2004). 

23. T. E. Hartman, "Tunneling of a wave packet" , J. Appl. Phys. 33, 3427 (1962). 

24. H. Winful, "Tunneling time, the Hartman effect, and superluminality: A proposed resolution of an old 
paradox", Phys. Rep. 436 1-69 (2006). 

25. S. De Leo and P. Rotelli, "Localized beams and dielectric barriers", J. Opt. A 10, 115001-5 (2008). 



Stefano De Leo et al: TRANSMISSION COEFFICIENT for quaternionic barriers 



9 




Fig. 1. Incoming energy dependence of the transmission probability for different values of the potential. The 
continuous line represents the standard case of a complex potential. The dotted line represents the case of a 
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Fig. 2. Barrier width dependence of the transmission probability for different values of the potential. The 
continuous line represents the standard case of a complex potential. The dotted line represents the case of a 



